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ERRATUM 


In the paper “The Mechanism of Cell Divi- 
sion” by N. Rashevsky, Bull. Math. Biophysics., 1, 
23, 1939, the following somewhat confusing nota- 
tion should be changed: substitute J for T in last 
term of equation (1); in left hand side of equa- 
tion (8) ; in first line from top on page 25; in 12th 
and 19th line from top on page 27; in left hand 
side of equation (20). 

In formula on 10th line from bottom on page 
24 the expression in braces should be multiplied 
by c. 


we 
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ELONGATION AND CONSTRICTION IN CELL DIVISION 


GALE YOUNG 
The University of Chicago 


_ The deformation of a constant reaction ellipsoidal cell by diffu- 
sion and constant surface tension forces is studied. The critical size 
of a spherical cell at which it becomes unstable to ellipsoidal defor- 
mations is found to be the same as that obtained previously by N. 
Rashevsky from energy considerations. It is shown that such a cell 
once unstable will elongate to a finite amount, and that it will tend 
to constrict in the center and round up at the poles. 


This paper makes an attempt to follow the deformation of a vis- 
cous body under the action of diffusion and surface tension forces. 
The treatment is in places admittedly rough and incomplete, but it is 
hoped that the work may present sufficient analogy with cell division 
phenomena to warrant its presentation. It is hoped also that it may 
serve to indicate the type of observations needed in connection with 
such theoretical studies. The dynamics of cell motion, which includes 
such phenomena as cell division, pseudopodic extensions and contrac- 
tions, protoplasmic streaming, etc., involves essentially the deforma- 
tion of a material system under the action of various forces. The sys- 
tem may offer both viscous and elastic resistance to the motion, while 
the driving forces may arise from a number of sources. Various writ- 
ers have assigned the driving force to various causes, usually however 
without trying to calculate the rate and nature of the motion which 
would ensue. It is clear that this is an essential part of any theory, 
and it is important to attempt to do this, even though the problem is 
mathematically and physically rather difficult. 

Among the more prominent views about the nature of the forces 
responsible for cell division is that of Butschli who regarded them as 
arising from a variation in surface tension over the cell surface, that 
of R. §. Lillie who considered them to be electrical in origin, and that 
of Heidenhahn who supposed that they arise from the pull of contrac- 
tile fibers (Wilson, 1934, pp. 189-198). It has since become clear from 
the work of N. Rashevsky (1938a) that the drag forces exerted by dif- 
fusing substances play an important role in the behaviour of meta- 


bolizing systems. 
27h 
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The present paper will consider the action of diffusion forces and 
constant surface tension. The first gives rise to volume forces through- 


out the system of amount 
F=—o Vc (1) 


where c is the concentration of the diffusing substance and a is a con- 
stant of the order of magnitude of RT/M, M being the molecular 
weight of the diffusing substance. The second gives rise to a normal 
pressure over the cell surface of amount : 


P=yK (2) 


where y is the surface tension and K is the mean curvature of the sur- 
face at the point in question. 

Instability. It has been shown by Rashevsky that this combina- 
tion of forces leads to a reasonable estimation of the size of living cells 
at which division might be expected to take place. First, by expand- 
ing the cell to infinite size and then recondensing it into two spheres 
each of half the original volume, it is possible to find a critical radius 
for the original cell above which the sequence of changes might occur 
spontaneously,:at least insofar as concerns the overall work done in 
the process. It is supposed that a varies inversely as the cell volume 
during the operation, and the same for the rate, q, at which the dif- 
fusing substance is produced per unit volume. For a very large mem- 
brane permeability the resulting value for the critical radius (Rash- 
evsky, 1938a, Chapter 8) is 


Ni 11.25 »y D 
oho be (3) 


ares ar 


nal and external diffusion coefficients D; and D.. 

Secondly, (Chapter 9) Rashevsky calculates the work done by 
these forces in infinitesimal deformations of the cell from the spher- 
ical shape. In particular we may consider, for reasons to appear 
presently, such a deformation which makes the cell into an ellipsoid 
of revolution. The equation in polar coordinates of an ellipsoid with 
semi-major axis a, semi-minor axis b , and focal separation 2f , is 


b e? cos? 6 
r= —________ = 0 il} See oe 
V1— e? cos? 6 ( an 2 sh ) (4) 


where e = f/a is the eccentricity. Making use of the relation 


2 
cos? 6 = 3 P,(cos 6) + P,(cos 6); 
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we find in Rashevsky’s notation that for e very small 
t.—0 
(5) 
1 
Ap a, = 3 b e?, 


while the other coefficients a, are negligible. Substituting these val- 
ues in his general expression (eq. 94, p. 97) and setting h = o, we 
find that the cell is unstable to such a deformation if its radius is 
larger than 


ce Dat SD) 


= (6) 


The difference between the two derivations should be carefully 
noted. Expression (3) gives an estimate for the radius at which the 
energy change from initial to final state becomes favorable to divi- 
sion; while (6) gives the radius at which the work balance is favor- 
able in a small deformation of the cell into an ellipsoid of revolution. 
They do not measure the same thing, but it may be noted that they 
agree well. If D = D, >> Dj,, then (6) gives a value 1.17 times 
that of (3); if D = D. = D;, then (6) is 1.39 times as large. 

Neither of the above developments concerns itself with interme- 
diate stages of division, or gives any assurance that division will pro- 
ceed. To consider these matters it is necessary to discuss in more de- 
tail the mechanics of cell deformation. The first stage of division in 
many cells is typically an elongation into an ellipsoid of revolution, 
the so-called “karyokinetic” elongation (Wilson, 1934, p. 157). This 
reaches a maximum value, and the cell is then divided by a furrow 
constricting inwards in its equatorial plane. Similarly the early stages 
of nuclear division are often marked by an approximately ellipsoidal 
shaped spindle, which later takes on a dumb-bell shape as the central 
region constricts. 

For consideration of the early stages of division it is thus appro- 
priate to study ellipsoidal shaped cells. This is why an ellipse was 
used in deriving (6). Considering the cell (or spindle) as an incom- 
pressible viscous body of viscosity 7 , and neglecting the resistance of 
the bounding membrane to the passage of the diffusing substance, it 
is found (Young, 1939) that the contribution to the average relative 
elongation along the major axis due to the internal diffusion drag is 


_ 2 aq . &—b? 
45 » 2Di+mD, ‘ 


€p 


(7) 
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where m is a function of a/b as shown on the accompanying graph. 
(Fig. 1). Similarly the contribution of the surface tension is 
y a—b 


€s.r. oie eee (8) 


where f is a function of a/b : 
BC) = 156 3 -6(@) = 153 BG) — Lis. 


Interpreting &@ = ép-++ és.r. just > 0 to mean the beginning of 
elongation, the critical radius is again found to be that given by (6). 
Thus there is complete agreement with Rashevsky’s energy treatment 
as regards the size of cell at which instability to infinitesimal ellip- 
soidal deformations sets in. With the present approach, however, it 
is possible to go farther and follow up to finite elongations — as far, 
in fact, as the cell approximates to ellipsoidal shape. 

Finite elongation. It will first be shown that if the radius exceeds 
the critical value given by (6), then the cell will undergo a finite 
elongation. In a constant volume ellipsoidal deformation 


ab? = constant = k = cae (9) 
An 
so that (8) may be written as 
An 
€s.r. =— 57, 6 ble). (10) 
Also (7) may be expressed as 
ate 0. b= 0/d 
TO EOS SARE OAT ome (11) 
Denoting 
2 D; SS ae Dy, “ (12) 
the ratio of (11) to the absolute value of (10) is 
Hed 1ta/b 1 
where 
ee V 
. 102yD.— (14) 


__. The first variable term in (18) has the value 5/4 at a = b and the 
value 2 at a = 2b. Instability occurs at r=1 3 1e., at 
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y == (e428). (15) 


To see if elongation will continue for this value of y it is necessary to 
compare 


1 b 
= (a % + 3) ae (16) 


etm. (17) 
As long as (16) exceeds (17) the elongation will continue. The for- 


mer has the value x + 3 ata = b, the value (a +3) ata = 2b. It 


and 


is thus approximately a straight line of slope S(t 3); which for 


comparison with (17) is to be drawn on the m-graph with the value 
at aN 
The minimum slope of the (16) line occurs for x = 0 (ie. 
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D> >-D;) and is 2 This is No. 1 on the graph and is seen to be 


above (17) up to a = 1.14 b, so that elongation once begun will con- 
tinue to this extent. 

As x increases the straight line swings upward, increasing the 
ratio of a/b at which elongation will cease. Up to position No. 2 this 
increase is gradual up to a= 1.3). Then, however, the equilibrium 
ratio suddenly jumps to about a = 5.76b as the hump in the m-curve 
is surmounted. It always remains finite, however, as will be seen 
later. This jump occurs at a slope of 2.76, which requires x = 1.6 or 
D; == 0.8 1D). . 

A line of slope between that of No. 8, which corresponds to 
x = 0.37 or D; = 0.185 D., and that of No. 2 will intersect the m- 
curve in 8 points. The middle of these is unstable, the top one stable 
(supposing, of course, that the formulae derived for an ellipsoid re- 
main valid). 

To see that the elongation (with the same supposition) would 
always remain finite, note that for a > > 6 the surface tension re- 


storing effect (10) varies as ab, i.e. as \/ka which becomes infinite as 
a increases. On the other hand (Young, 1939, p. 49) as the elongation 
2 


; 2 ‘ Oo re 
continues €p varies as b? log — » Le. as 
a 


6 


b? log g 


ke (18) 


which vanishes as b tends to zero. 

Thus as far as ellipsoidal shapes are involved, a cell once unstable 
will continue to elongate to a finite amount. One concludes that the 
oe is probably true for other shapes not too far removed from ellip- 
soids. 

A somewhat similar conclusion is reached by considering instead 
of (7) and (8) the corresponding formulae obtained by Rashevsky by 
an approximate method for a more general class of shapes (Rashey- 
sky, 1939). If the quantity (a + 1) /2 is greater than unity, where 
x is as in (12), his expressions give the result that if é is just greater 
than zero at a just greater than b then it vanishes ata = («7 -+1)b/2. 
Between these two values for a/b the quantity é is positive, so that 
when instability to small elongations sets in the elongation will con- 
tinue to a finite amount determined by the ratio of the diffusion co- 
efficients. Thus in this case the conclusion agrees with that obtained 
above from the study of ellipsoidal cells. For x = 1.6 both treatments 
give a = 1.3 b for the amount of elongation. 
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On the other hand if « < 1, i. if Di; < 4 D., then the other 
equilibrium ratio in Rashevsky’s treatment, namely a = (+ 1)0/2, 
corresponds to a flattened or shortened shape body. When the spher- 
ical shape is just unstable to small shortenings the shortening will 
continue to a finite amount. The above treatment of ellipsoids has been 
confined to prolate or elongated ellipsoids, and hence can not be used 
to discuss the case of b > a. Extension to oblate or shortened ellip- 
soids would permit comparison with Rashevsky’s results in this case 
also. This will not be entered into here. 

Rate of elongation. If D; < 0.8 D., which will usually be the 
case for entire cells though perhaps not for mitotic spindles, then a 
reaction rate just large enough to make the spherical shape unstable 
will not carry the ellipsoidal elongation beyond a = 1.3 6b at the most. 
In some instances, particularly in division of entire cells, the early 
elongation does remain within some such small value. In other cases 
however, and usually so with spindles, the elongation exceeds this 
value while still remaining approximately ellipsoidal. If D; < 0.8 D. 
and the elongation does progress beyond a = 1.3 D, then the reaction 
rate, if supposed constant throughout the cell, must have increased 
since the beginning of the elongation.* 

In any case the average rate of relative elongation at any mo- 
ment is proportional to é = ép + és.7,, i.e. to the height of the straight 
line above the curve in Fig. 1. If the reaction rate is considerably 
larger than that necessary to instability, this height will be consider- 
able and as an approximation we may neglect the bends in the curve. 
Over a considerable range we may then replace m by 3a/b in (7). 
This quantity is somewhat greater than 2, and if D, is considerably 
greater than D; we may neglect the latter in the denominator of (7). 
Similarly we shall ignore the small change of fin (8) and set 6 = 1.5. 
Identifying é with (1/a) (da/dt) we then have 


y a—b 
ee a 10?) 9 19 
a dt 135 Dea ‘ ) 2 ab (12) 


It may be that in a given set of data the variation in b? is small 
compared to that of a? , and one might replace b* by ka? in the a? — b? 
term, where k is an average of the ratios at the ends of the data range. 
For example, in a set of data to be considered presently b? varies be- 


*In the present paper we are supposing the reaction rate constant through- 
out the cell, but it is not insisted that this is the case. There are in fact reasons 
for thinking that astral systems represent mechanically important localized re- 
action centers, and one would surely want to study them as such. Here, however, 
“we are tracing consequences of the constant reaction assumption. 
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tween .4 a? and .3 @?, and we shall take k = 0.35. With this simpli- 
fication (19) reduces to 


fe RAR SOLNC Seni 8; wo (p—1). (20) 
dt 185 n De 2n 
In a constant volume ellipsoidal deformation 
an a b? = V = constant , (21) 
and (20) becomes 
da 4) 
re ee >, De 
TE =A—Ba (22) 
where 
reich te 
and 


ee (Aaa PnP ed gilt (24) 
2n'3V 185 7D. 4n 

We shall attempt to apply these results to some observations of 
Béla¥ (1929) on the elongation of mitotic spindles. There are a num- 
ber of objections to this application, of which but two will be men- 
tioned. The first is that the spindle volume is not constant in Bélar’s 
material, so that the increase in length comes about from growth as 
well as stretching, which would tend to exaggerate the velocities. The 
second is that the spindle is not a homogenous viscous body. Follow- 
ing Rashevsky (1938b) we shall attempt to reduce this difficulty by 
considering only his data after the chromosomes have separated, 
which should decrease to some extent the effect of possible contractile 
or elastic fiber pull. 

Bélar’s data gives spindle length as a function of time. It was 
first smoothed by fitting a curve to it (Rashevsky, 1938b, Fig. 5), 
and then the values of a and da/dt were taken from the curve. For 
any two points (22) gives by subtraction 


so that each pair of points gives a value for B. For simplicity each 
other point was paired with the second point, and the resulting values 
averaged to give B. Using this average B in (22) gives an A value 
for each point, and these were similarly averaged. The results are 
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given in the accompanying table, where the last column is obtained 


by substituting the average B and A values back into (22). These 
values are 


A=1.177X 10" 


(26) 
Boe 0.050. 
observed 2a smoothed 2a | smoothed da/dt |calculated da/dt 
in 10-+ cm. in 10-4 cm. |in 10-6 cm/sec Jin 10-6 cm/sec 
21.4 21.8 0.51 0.48 
24.4 24.3 0.22 0.24 
VAS 25.4 0.11 0.15 
25.6 25.9 0.083 0.10 
26.6 26.4 0.07 0.04 
27.0 27.0 0.008 —0.03 


The volume V is obtained by averaging the initial and final val- 
ues, much as was done above for k ; 


V =2.6 10° 
(27) 
Be 0Sh: 
Using these values in (23) and (24) we obtain 
Y = 3.54 10 (28) 
| 
and 
7 14105. (29) 
n De 
Taking 
Lins 
com 
M = 100 (30) 
Dee 1054 
gives 
4 = 5.7 10 (31) 
| 


Using (28) and (29) in (6) with D. > > Di gives for the criti- 
cal radius 7.7 < 10+ cm., while the volume of (27) corresponds to 
a radius of 8.57 « 10 cm. for the spindle in spherical shape. As 
mentioned above this would, in the present treatment, indicate that 
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the reaction rate is greater than required to cause instability of the 
spindle in spherical shape.* 

Constriction. The above considerations hold only insofar as the 
cell remains approximately ellipsoidal in shape, and their application 
depends upon the empirical fact that cells and spindles do have ellip- 
soidal stages. So far the theory has accepted the shape as given, 
without trying to see what sequence of shapes would actually be ex- 
pected. By the use of Betti’s formula (Young, 1939) the average rate 
of elongation of an incompressible viscous body under the action of 
any assigned forces can be determined, as is illustrated by expressions 
(17) and (8). This however gives little information as to the change 
in shape of the body; it might be elongating in some regions and con- 
stricting in others. For application such as to cell division it is nec- 
essary to study different regions of the body in greater detail. 

To this end consider a slab of the body cut out by two parallel 
planes. If we can estimate the average rate of elongation in such an 
arbitrary slab we shall have more information about what is happen- 
ing in different parts of. the body. Again considering only the effect 
of internal diffusion drag and constant surface tension, the forces act- 
ing on the shaded portion of the body (Fig. 2) are (a) surface ten- 
sion pressure over its lateral surface S; ; (b) its internal diffusion 
drag forces; and (c) the stresses exerted across S, and S, by the end 
regions of the body. 


FIGURE 2 


In determining by Betti’s formula the average elongation é of the 
slab in the z direction the effect of the drag forces (b) is the same 
(Young, 1939) as if a pressure ae were acting normally over S, + 
S.-+ S; Assuming for simplicity that the system is symmetric about 
the z axis, the (c) forces involve the integral . 


SS (2Z’—a«X’')ds, (32) 


81482 


where Z’ and X’ are the stress components exerted by the end regions 


* This of course does not imply th i 
; at th i i 
shape, it may have come into beitig ti an slatoutea Hae nan ‘ ateca, Biediiaeaa 
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across unit area. The Z’ part is readily evaluated since z is constant 
on S, and on S,. Thus 


SfzeZds=2,ffZds=2,F, 
81 


(33) 
where F, is the total z force acting on the right hand end region; 
namely F’, = (total diffusion force on the end) — (lateral surface 


tension holding the end on). The first term of F, is as if a pressure 
ac were acting over S, + S,. The S, part of this cancels the S, part 
of the (b) term, and similarly for S,. Thus we are left with the fol- 
lowing contributions: 


(I) pressure a ¢ acting on S; 


(II) 2 [(z force of pressure a ¢ on S,) + 
(z surface tension pull on end no. 1) ] 


(III) 2. [ (2 force of pressure ac on S;) + 
(z surface tension pull on end no. 2) ] 
(IV) —ffuX'ds 


81482 


(V) surface tension pressure acting over S;. 


Note that a constant added to ¢ will cancel out from (I), (II), 
and (III), because in (II) a constant normal pressure over S, gives 
the same total z force on end no. 1 as if it were acting normally (from 
the right) over S,, and similarly for S; and S,. Then with (I) the 
effect on the slab is of a constant normal pressure over its entire sur- 
face, which contributes nothing to @. Thus only variable terms need 
be retained inc. 

Constriction in ellipsoid. We shall first apply these considerations 
to an ellipsoid of revolution to see what is the tendency for departure 
from that shape. In this we shall make use of notation, methods, and 
results previously given (Young, 1939). The longer axis is along z ; 
the semi-major and semi-minor axes are denoted respectively by a and 
b ; the distance of the foci from the center is f = ~@ — 0 sand 52 
_ distance from the center is z = a cos 6. On the ellipsoidal surface the 
concentration is of the form 


c=G,+ G cos? 6, (34) 
wherein, as mentioned above, the G, term may be ignored. The con- 
tribution of term (1) to é@ is 


a 

ego ; | (35) 

where 

J = Sfcl[z(cos nz) —2xcos(nx)] ds, (36) 
Bs 
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and where V is the volume of the slab under consideration. Carrying 
out the integration gives 


J=n2Ga b (SL cos" a:—* [cos® 6]?} 


(37) 
7 b? a am b? 
=56— (2%, — 2°,) Steger A keh (25, — 21); 
where the subscripts 1 and 2 refer to the two ends of the slab. 
The total z force of pressure ac acting over S, is given by 
—affecos(nz) ds; (38) 
S4 
which, apart from the G, term, turns out to be 
5 Ges eee nee. (39) 
2 Jewey) (ee 


The z component of surface tension holding end no. 1 onto the slab is 
in2 
phaendt ab sin? 6, (40) 
Va? — f? cos?6, 


Similarly the total z force acting on end no. 2 is the sum of 


b2 
506 (a — z4,) (41) 
and 


a? — f? cos? 0, 


Finally the contribution to é of surface tension over S; is found 
to be die 


ahr Sao ie aS (43) 
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3 Lae eae 
35 pe AV — f? 27; — 2, at — fF? 25) ]- 


These expressions simplify considerably if the slab is made very 
thin, so that z, — z. = dz. Thus (37) becomes, by differentiation, 
ae 2 
= 2 BE pas Ti a ° 
J=21Gb0?(3 2 Ge) a 5 (44) 
which with (85) and 


Sp ei ton (45) 


gives a contribution 
CS ae ee ee (46) 


The terms (39) and (41) contribute 


y 1 
e= —_. [z, (89).-— 22(41) 1. (47) 
SPV. 
The value of the bracket comes out to be 
n b? 
oe Ge [2°; — 2 at + 2 at — 2°], (48) 


which for a thin slab reduces to 

n b? 

— ——— 4 —_—__ A Key @ 

gt Ge (52 a‘) dz (49) 
In this last case (47) becomes 


° Bn 2a a?— 2 ee) 

The term (IV) can not be directly evaluated from a knowledge 

of the forces applied to the entire body, as the X’ relate to the state of 

internal stress set up. For the present this term will be ignored. With 

this neglected the total é in the slab due to diffusion forces is the sum 
of (46) and 50), namely 


CN cy ee (51) 


Similarly (43) reduces, by differentiation, for a narrow slab to 
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1 
eS Oe. — Ot 
é = ab dz ara fray | + 
(52) 
(6 at — 2 a? f?) 2? — 8 fF 24] ; 
while the term 
1 
oc 1 2 42 53 
é 3,v (40) + 2, (42) ] (53) 
becomes 
= — 2b ae —payn(a—Bare+2fe). (54) 
| 


Adding (47) and (49) and introducing (45) gives for the total sur- 
face tension term, (IV) being neglected as before, 


ég.r.(2) = Sy he Peyr: (55) 

It is seen that both (51) and (55) are numerically largest at the 
equator (z = 0) and drop off to zero at the poles (2 = + a). They de- 
crease, however, at different rates, and this gives rise to an interest- 
ing consequence. With G negative expression (51) is positive, show- 
ing a tendency for each part of the body to elongate; while (55) is al- 
ways negative, showing that surface tension tends to thicken each 
part of the body. But (51) drops off toward the poles more rapidly 
than does (55), so that the thickening or rounding up tendency is 
relatively greater at the poles, while the elongating tendency is rela- 
tively greater at the equator. Since the neglected term (IV) involves 
transverse stresses which are opposite on two adjacent slabs, its re- 
tention would tend to smooth out the variation in €(z) from poles to 
equator. 

In particular if é@ for the entire body is zero, then there will be 
constriction at the center and thickening at the ends. Expression (7) 
is (Young, 1939) 


‘eae 2aG 
De mebon 4 ? (56) 
and thus (7) plus (8) equals zero if 
_ d5yBba—b 
Ca: Oe Sanne re 


With this value of G the sum of (51) and (55) is 


(a? — 2?) —" a? (at b) 5B 


Mh Taortae nwa woraie| 8) 
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which changes sign at cos 6 = z2/a satisfying 
4 
5p 


where e = f/a. The central region of the cell up to values of z given 
by (59) will be constricting and elongating; the end regions for great- 
er values of z will be broadening and shortening. 

For example consider the case a = 2b and thus e? = #,8 = 1.5. 
Then (59) gives z = + 0.43 a, so that the ellipsoid is constricting 
over a little less than half its length, and rounding up over the re- 
maining. If a = 1.34 b, with e? = 4/9 and 6 = 1.57, then zg = + 
0.31la. 


e? cos? 6 = teat, a+2] , (59) 


: 1 da : : 
In the above, as in (19), snap was identified with the volume 


average é over the entire cell interior. Instead of employing (7) or 


(56) as the diffusion contribution to - i , one might use 


5 oe pe 

re =< | ép(z) dz, (60) 
a Jo 

with é@(z) as given by (51). This comes out to be 

ép = —=— > (61) 


which differs slightly from (7) in the numerical factor. In the same 
way one might consider, instead of (8), 


és a he aT: a 1 f 
ae = —_—_ — ~~ : 2, 
an = 5) eax. (2) dz a(t asin ) (62) 


This may be written as 


iy ae (68) 
o 7 ab 


where £ is a function of a/b ; € (1) = 1.33, € (2) = 14,& (0) = 


1.57. Using (62) plus (61) gives the same critical radius (6) as using 
(7) plus (8), and in general the results are closely comparable. 
Constriction in other shape bodies. It has been seen above that 
as the body elongates into ellipsoidal shape it tends to pinch in at the 
equator and bulge at the poles, which is reminiscent of the actual 
course of cell and spindle division. As soon as the constriction has 
proceeded to any appreciable degree the above formulae cease to hold, 
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and must be revised. It does not appear feasible to attempt to follow 
the complete course of the deformation, but we may make some rough 
estimates. Suppose the constriction to have proceeded to the stage 
of the diagram (Fig. 3), where the body is comprised of two spheres 


FIGURE 3 


connected by a narrow neck. The diffusion problem is now somewhat 
like that of two separate spheres, and thus (Young, 1936) the force 
tending to separate the spherical ends is roughly given by 


fee 2s Oa eR (64) 


The lateral surface tension force tending to hold each end onto the 
neck is of the order of 


ZITA os (65) 


while over the lateral surface of the neck itself is a normal surface 
tension pressure of about 


? 
r (66) 


As the spheres tend to pull apart the narrowest region of the neck 
constricts inward, much like the “necking” of an iron rod in a testing 
machine. Introducing the above forces into Betti’s equation for the 
central part of the neck alone, and ignoring its internal diffusion 
forces, we obtain 


1 F—ary 


SS Say 7 (67) 
Now 
Pio PMO 
"medial ia 


so that if a and b be regarded as remaining roughly constant while 
the constriction proceeds we have 
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where 
M = y/6n 
(70) 
| bé 
36 7 D, a 
At this rate r would drop to zero in a finite time, provided that it has 
started to decrease at all. If N/r > > M, the decrease is according to 


dr N 

tena (71) 
so that from a finite initial value 7) the radius will have dropped to 
zero in a time. 


ine 
2N- 


t= (72) 
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ON THE SO CALLED PLASTICITY OF THE CENTRAL 
NERVOUS SYSTEM 
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The taking over of a function of a region of the brain by an- 
other region, in case of brain injuries, is discussed in its mathe- 
matico-biophysical aspects. 


One of the most important properties of the central nervous sys- 
tem is a peculiar combination of definitely localized structures on one 
hand and of a large amount of variability of function of a given struc- 
ture on the other. Removal of a region of the brain, which is known 
to control a given activity, usually results in the loss of that activity. 
But not infrequently the loss is temporary, the lost function being 
eventually assumed by some other part of the brain, which normally 
may not appear to contribute anything to that particular function. It 
is the purpose of this note to point out, that such phenomena find a 
very natural and simple explanation in the light of the mathematico- 
biophysical theory of the brain, developed by the author in his book 
(Rashevsky, 1938a) and by the author and others in subsequent pub- 
lications (Landahl, 1938, Rashevsky, 1938b, Householder, 1939). 

In these publications the following general structure has been 
considered and applied to various special cases. A number of parallel 
excitatory neuronic chains are inter-connected with each other by in- 
hibitory fibers, so that the excitation of one neuronic chain tends to 
inhibit the others. As has been shown (Rashevsky, 1938 a, b, House- 
holder, 1939) such an arrangement possesses the following property: 
If all the parallel chains are excited with unequal intensity at their af- 
ferent ends, then the excitation is transmitted only through a few of 
these chains, which are excited most strongly, all the weaker ones be- 
ing completely inhibited, even though they are stimulated with a 
super-liminal inbensity. How many of the chains will be inhibited 
depends on the relative values of excitation of the different chains, 
and on the constants characteristic of the inhibitory fibers. 

Suppose now that a given end organ is innervated by a number 
of central neurons, spontaneously excited (Rashevsky, 1938a, Chapt. 
XXIV), which are scattered throughout the whole cortex, or through a 
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very large region K, of the latter. Let however these neurons be not 
identical in their physical properties, but possess different intensities 
of excitation. Let the most strongly excited neurons lie more densely 
around a point P(ao, Yo, 20), which is contained in Ky. The density 
of the more strongly excited neurons may fall off continuously from 
a certain point P to the other parts of K,. We may say then that to 
each point x, y, z in K, corresponds a given average intensity of total 
excitation E= f(a,y,z), where f(x,y,z) has a maximum for % = %, 
Y = Yo,% = %. If each of these neurons sends off inhibitory fibers to 
all others, then under these conditions (Rashevsky, 1938b) all neu- 
rones, having an intensity of excitation EH less than a certain value E,, 
will be inhibited, while those with HF > E, will remain excited. But 
those latters lie in a sharply limited region K, around P, the bound- 
ary of this region being given by f(2,y,z) = H,. Application of a 
stimulus to region K,, which will change the excitatory state of the 
neurones contained therein, will result in a change of the state of the 
end organ. The application of a similar stimulus to any point outside 
of K,, but within K,, will not result in any change of state of the end 
organ, because of the complete inhibition of the neurones involved, 
by the more strongly excited neurones of K,. Thus, such a “func- 
tional” test will reveal a sharp localization of function. Let however 
the region K, be removed. The inhibitory action of the neurones con- 
tained in K, upon the remaining ones is thus also removed and some 
of these remaining neurones will become now excited and appear 
functional. Again however not all of the remaining neurones will 
carry on the excitation, but only a small group of the most strongly 
excited ones. This group is again localized in a sharply defined region 
K,’. If the function f(a,y,z) falls off monotonically with increasing 
distance from P , then the functions of a removed region will be taken 
over by an adjacent region. But if f(x,y,z) has several relative maxi- 
ma, K,’ may lie quite far apart from K,. In general we should expect, 
that the extent of damage to K, would affect the size of K,’' , depend- 
ing on the shape of the function f(a,y,z). It should be possible in 
principle to determine the shape of f (x,y,z) by a quantitative study of 
the effects of brain injuries of different extent. 
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_ Various neural mechanisms are considered which deal with 
point to point correspondence between two sets of neural elements 
with a smaller number of conducting elements between them; the 
transmission of nerve impulses in a limited range of intensities; 
movement of the transmission of excitation along a contour; the 
reaction to the size of an object independent of its distance; and an 
interpretation of the effect of a warning stimulus and of stimulus 
intensity upon reaction time. For the latter cases a comparison of 
the theoretical equations is made with some of the available experi- 
mental data, and a general agreement is found. 


In discussing the following nerve mechanisms we shall use the 
concepts and notation developed by N. Rashevsky (1938, Ch. XXII), 
unless specific reference is made to the contrary. If the afferent proc- 
esses of a nerve fiber are stimulated by an intensity, S, greater 
than a threshold, h., excitation occurs and the intensity of the excita- 
tion, E (proportional to the frequency and strength of the individual 
nerve impulses), increases in some manner with the strength S of the 
stimulus. At the efferent processes of the nerve fiber, an excitatory 
factor, «, or inhibitory factor, 7, (or both) is produced at a rate pro- 
portional to the excitation, EF, and is dissipated proportional to its 
own concentration. The difference « — j acts as a stimulus on the 
next order nerve fiber if e — j is greater than the threshold, h, of the 
latter fiber. 


I. Point to point correspondence between two sets of neural ele- 
ments with a smaller number of conducting elements between them. 
It is well known that there are more photo-sensitive elements on the 
retina than there are transmitting elements in the optic nerve (Maxi- 
mow and Bloom, 1938, Ch. XXVII). Only certain, cones are con- 
nected through bipolar neurons to ganglion cells in a one-to-one rela- 
tion (b—h—s), though the third order neurons are influenced by other 
cones (Maximow and Bloom, 1938, p. 601). If a group of rods or 
cones stimulates a single ganglion cell and the result is only an in- 
creased sensitivity, no difficulty would arise. However, some evidence 
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indicates that a higher degree of discrimination is possible between 
elements in the retina than would be indicated by the number of fibers 
in the optic nerve. The net result would be that a group of neurons 
is able to transmit excitation to corresponding members of another 
group through a single fiber or smaller number of fibers than in the 
initial or final group. Such a mechanism could conceivably occur else- 
where than in the optic system. For convenience we shall discuss the 
problem as though the phenomenon occurred in the optic system. 
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FIGURE 1 


An interpretation of such a phenomenon must not involve the 
use of hypotheses which are contradicted by known facts regarding 
neurons and their behavior, or the use of specific assumptions, unless 
the latter are indicated directly by more fundamental processes. One 
might otherwise assume a specific type of action by each member of 
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the first group which could be transmitted only by the corresponding 
member of a higher group; the specific type of action might be a cer- 
tain frequency, intensity or pattern of both. The well known work of 
P. Weiss (1936) seems to indicate specificities in groups of fibers. We 
shall, however, attempt a solution without the use of such an assump- 
tion. 

Let us now discuss the mechanism in Figure 1 which exhibits the 
property of apparent point to point correspondence between two sets 
of nervous elements with a smaller number of intermediate connect- 
ing fibers, but which involves neurons having properties assumed from 
other lines of evidence. The mechanism was suggested while study- 
ing Dr. S. Polyak’s semi-schematic diagram of the retina (Maximow 
and Bloom, 1938, Fig. 519), although there may seem to be little in 
common between Figure 1 and Dr. Polyak’s diagram. 

In order to have a figure which is not entirely too complex for 
illustration, we have in Figure 1, five groups of only four fibers (A1, 
A2,---,B1,---). Each of the four fibers in a group make synaptic 
contact with a common fiber (A , B, C---) which in turn makes synap- 
tic contact with each of four third order fibers (Al’, A2’,---, Bl’, 
.--). Let each first order fiber (A3) also make synaptic contact with 
an additional afferent fiber (B) which invariably is a collateral from 
a different group, and let this collateral afferent fiber have a sym- 
metric collateral efferent fiber which connects with the third order 
fiber (A3’) corresponding to the first order fiber. Let us also assume 
that the threshold of the third order fibers is such that excitation 
does not occur unless both efferent fibers which synapse with it are 
excited together. For convenience we assume that all fibers, having 
the same function, are identical, though it is not necessary (Hecht, 
1927-1928). Also the omission of the bipolar cells from consideration 
may have introduced an. unnecessary restriction. 

It is possible to satisfy the required conditions in at least two 
ways. In Figure 1 a closed system is shown, and it is also possible to 
construct similar systems which are complete to any multiple of the 
above. That is, with four fibers in a group, a complete system would 
involve forty first order fibers, ten second order and forty third order 
fibers, the ratio still remaining four to one. Let us consider only the 
simple case shown in the figure. The efferent end of every first order 
fiber, as A3, connects with a second order fiber of its own group, A, 
and an afferent collateral from another group, B ; and no two fibers 
in the same first order group synapse with a collateral from the same 
second order group. The cyclic arrangement in the figure is quite evi- 
dent upon inspection. If we refer to A as being above B, then we 
may state that the fourth fiber of the first or third order connects with 
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the second order fiber of the second group below, as A4 —C , B4—D, 
C4 — E, DA — A, E4 — B, with similar connections between the 
second and third order, C— A4’, D — B4’, etc. In a similar manner 
the third fiber of the first or third order connects with the second or- 
der of the group directly below, as A8 — B, B38 —C, C3 — D, ete. 
And if we take a group of first order fibers, as the C group, we have 
Ci —A,C2—8B,C3—D,andC4—E. 

Under the above conditions, a fiber of the first order, as A3, can 
then excite its corresponding third order fiber, A3’. However, the lat- 
ter, A3’, may be excited by the simultaneous stimulation of any fiber 
of group A together with any fiber of group B. In order to meet this 
objection we shall consider the mechanism from two points of view. 
First we consider all fibers identical and also that the time during 
which two impulses can summate at a synapse, or the period of fa- 
cilitation, is very short (Lorento de N6, 1938). Then we shall assume 
that there are two types of fibers and also that properties of the fibers 
are as discussed by Rashevsky (1988). 

If the period of facilitation at the synapses between the second 
and third orders, as A — A3’, is small, then the two impulses initiated 
by the first order fiber, A3 , arrive together at the third order synapse, 
A3’ (we assume that the second order fibers are very long as com- 
pared with any collaterals), and by summation can exceed the thresh- 
old of A3’ (Lorento de N6, 1938) ; whereas, an impulse initiated by 
Al, A2, or A4 and an impulse initiated independently by a fiber in 
the first order of group B would not arrive at A3’ within the summa- 
tion period, except under rather special circumstances, and hence the 
summation would be insufficient to exceed the threshold of A3’. Ina 
similar manner B2’ would be excited by A3. If either A3 or B2 are 
excited both A3’ and B2’ will be excited. If, in the first and in the 
third order, two fibers constituting a pair (as A8, B2) are together, 
a group with such pairs taken as units will show the properties re- 
quired — a correspondence between units in the first and third orders: 
a separate study would be required to investigate the special conditions 
under which independently initiated impulses would summate, and 
more specific assumptions regarding the mechanism and the actions 
of the various fibers would be required. 

Let us now assume that the period of summation is long compared 
with the period between impulses even at low intensities of stimula- 
tion. Then let the afferent and efferent ends of the second order fibers 
be of two kinds. Let those which connect with another group (light 
lines in Figure 1) be capable of carrying a greater excitation (i.e. hav- 
ing a higher limiting value of E' ; Rashevsky, 1938, p. 218) than those 
which connect with the same group (heavy lines). In fact, assume 
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that the simultaneous maximal stimulation of all of group B together 
with A4 will not quite excite A3’ , whereas the stimulation of A3 alone 
will produce a sufficiently strong excitation in fiber B which is able to 
summate with a weaker excitation through A at the synapse to exceed 
the threshold of A3’. Under these conditions there will be a point to 
point correspondence between the first and third order. 

One-half of the mechanism of Figure 1 is drawn in Figure 2 at 
right angles to the former figure and slightly rearranged. The con- 
nections are the same as in Feure 1 as is also the labeling. The solid 
circles represent the second order fibers whereas the open circles may 
represent either the first or third order fibers. The figure is then in 
the plane of the retina or “central retina” (Rashevsky, 1938, p. 266) 
and only one complete unit is shown. Though the connections seem 
rather random at first glance, the same cyclic symmetry is evident 
when it is looked for. 

Whether or not the above principle is found to exist anywhere 
in the central nervous system, it is worth while to note that a very 
simple symmetric pattern, as in the figures, may be very obscure and 
would be difficult to observe unless it was looked for, after being in- 
dicated by some hypothetical consideration. 

Summary. A particular symmetric arrangement of nerve fibers, 
which have properties customarily attributed to them, is shown to ex- 
hibit the property of point to point correspondence between two sets 
of excitable elements when the number of intermediate fibers is small- 
er than the number in either of the connected sets. 


II. A mechanism transmitting nerve impulses in a limited inten- 
sity range. Let a stimulus S be applied to two nerve fibers which 
arise and terminate together (I of Figure 3) ; one excitatory (solid) 


I 
FIGURE 3 


and one inhibitory (dotted). Let the thresholds for the two fibers be 
equal and equal to h. The intensity of excitation in the excitatory 
fiber, H., and in the inhibitory branch, Z;, may be given by (Rash- 
evsky, 1938, Ch. XXIT) 


hs =—. [kee 90(S-h) } (1) 


and 
ee oe [1 eas 8] (2) 
j 


100 MATHEMATICAL BIOPHYSICS 


where J is the intensity of each individual nerve impulse, @ is the re- 
fractory period and « is a constant of proportionality. 

At the synapse between the pair of fibers and the next higher or- 
der neuron, s, the excitatory and inhibitory fibers produce respec- 
tively the excitation factor ¢ and inhibition factor 7 according to the 
equations (Rashevsky, 1938, Ch. XXII) 


de 
= = AE —e (3) 
dj 
di =BE;—)j. (4) 


Solving equations (3) and (4) and imposing the initial conditions 
that e = 7 = 0 for t = 0, we have 


€ aes (1 —e*') (5) 
ets ey (6) 


In order for excitation to occur in the higher order nerve fiber II, the 
net excitation must exceed the threshold, e — j > h’, or using equa- 
tions (1), (2), (5) and (6), we have, 


We< A eas ae —— @-4e mic (<= es) 


(7) 
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Of several possibilities, let us here consider the case that 
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a eg 77 b 6; 


These conditions may be written in a somewhat more special form: 


Aclt. JeB i: 
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Under conditions (8), relation (7) may be written with close approxi- 
mation for the stationary state 
Al, a8 ais 


eh: a 6. 


~(S—h). (9) 


Hence, in order that the second order fiber be excited, we must have 


H. D. LANDAHL 101 


the stimulus intensity, S , lie within the range 


b AME ; 

Ti dd (10) 

which can be satisfied if 
mle 11 
or (11) 

Let 
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so that equation (10) may be written 
ht+é6>S>h (18) 


and only stimuli whose magnitudes lie between h and h + 6 are ca- 
pable of eliciting excitation in the second order fiber II and that, only 
if 6 is positive. A similar result as this was obtained by postulating 
a more complex mechanism by N. Rashevsky (1938a). 

Now if 6 is very small, values of S approximately equal to h, are 
the only ones which produce excitation in II. The width of the range 
of intensities transmitted, 6, decreases with increasing h’. This is 
equivalent to stating that the acuity of discrimination increases with 
increased inhibition (increased threshold) at the synapse between I 
and II. 

We have treated the fiber pair I in Figure 3 as two fibers, one 
excitatory and the other inhibitory. Evidently they could as well be 
a single fiber producing two factors, e and j , which are not only pro- 
duced at different rates but which also depend in a different way upon 
the stimulating influence. The results which have been obtained from 
conditions (8), in particular, would be unchanged, so that such a gen- 
eral type fiber may, under certain conditions, transmit excitation with- 
in a limited range of stimulus intensities. 

Equations (1) and (2) are of a specific form, but evidently more 
general relations could have been used. If we designate these rela- 
tions by E.(S) and E;(S), then, in order that transmission will occur 
in a single, finite range we must have [equations (5) and (6) for 
ae | 


Boy aig caihenelae MER (14) 
a b 


in a continuous finite range, h abe 6, > S >h-+ 6., where h + 6, and 
h + 6, are the solutions of (14) with the inequality sign replaced by 
the equality sign. 
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Summary. If a fiber producing only the excitation factor, ¢, ex- 
cites a fiber of the general type with the properties discussed above, 
then for the stationary state, only a certain range of frequencies in 
the first order fiber will be transmitted through the general type fiber 
to some third order fiber, since ¢, of the first fiber is determined by 
the excitation E, of the first fiber, which is by assumption (Rashev- 
sky, 1938, Ch. XXII) proportional to the frequency of the impulses. 
Another first order fiber, having different parameters, would have a 
different range of frequencies transmitted. And if the two fibers act 
together, only certain combinations of frequencies could be trans- 
mitted. 


III. A mechanism in which transmission moves along a contour. 
Consider now a large number of identical fiber pairs whose afferent . 
ends are excited by the efferent ends of one or more lower order neu- 
rons. Let these points of contact occur in a random manner so that 
there is a distribution in the number of synapses between the lower 
order and higher order fibers; and, as a result, the higher order fiber 
pairs will be excited variously by intensities ranging from zero to 
some maximum value. Then only those fibers whose stimulus inten- 
sities range from h to h + 6 [equation (13) ] will transmit the excita- 
tion. If now the initial stimulus applied to the lower order neuron is 
increased, all the intensities exciting the fiber pairs are increased and 
a different set of fibers will transmit the excitation due to the inten- 
sities ranging from h toh-+ 6. 

On the other hand, we may consider a large number of similar 
fiber pairs, differing only in respect to the threshold. Then, if they 
are equally excited by some stimulus, only those will transmit excita- 
tion whose synaptic thresholds range from h — 6 toh ; or from h, to 
ho + 6 if ho =h—o6. The net result here is the same as that obtained 
by the mechanism discussed in the preceding paragraph. More gen- 
erally, then, we may have a variation in both the number of synaptic 
contacts between the lower and higher order fibers and the threshold 
values of the fiber pairs. 

If the afferent processes of a large number of fiber pairs (Type 
I) of various thresholds are in synaptic contact with the innumerable 
collateral branches of the axon ending of the lower order neuron (er. 
neuropil), we have the more general case discussed above. Let the 
lower order neuron be of the intermediate excitatory type (Rashev- 
sky, 1938, p. 260) and let its excitation and inhibition constants a’ 
and b’ be very much smaller than either a or b of the fiber pairs [equ- 
ations (3) and (4)]. The fiber pairs are then always ina nearly stable 
position as the net excitation, (e—j), in the synapses between the 
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first order neuron and second order fibers rises and then falls. We 
may then use equation (13) with S replaced by « — 7 since e — 7 plays 
the role of a stimulus. But « — 7 first increases from zero rapidly to 
a maximum, then decreases somewhat more slowly to zero, so that 
the point of maximum net excitation (« — 7) is nearer the initial time 
than the second time at which e—j7 = h,. 

If we consider the fiber pairs arranged in the order in which ex- 
citation occurs with increasing stimulus strength (for one to one 
synaptic connections, in order of increasing thresholds h ; for equal 
thresholds, in order of decreasing frequency of synaptic contacts), it 
follows from the previous discussion that, when the first order neu- 
ron is given a constant stimulus, at first the low threshold fiber pairs 
transmit the excitation for a short time and subsequently higher 
threshold fiber pairs transmit. The excitation passes over as a wave 
of variable width from the fibers of lowest thresholds to the fibers 
whose thresholds are (€ —7) maz; the wave of excitation then passes 
back to the lowest threshold and stops, the return being slower than 
the forward process. If the secondary fibers are arranged at random 
in space, the appearance of a wave as such would be entirely obscured. 
The result would be more analogous to a great many lights, each flash- 
ing on momentarily in no apparent order and only slightly overlap- 
ping each other in time. The last few to come on would remain on 
longer; then a complete reversal would take place — the lights flash- 
ing in reverse order and at a slower rate, with longer duration of the 
light periods. 

An interesting case arises if the secondary fibers are arranged in 
order along a contour and if the first order neuron exhibits a slow 
periodic fluctuation of the excitation factor. (Landahl and House- 
holder, 1939). The transmission of excitation then travels as a wave 
back and forth along the contour with possible delays at each end. A 
mechanism with this property was suggested by N. Rashevsky 
(1938a), in connection with the development of a theory of aesthetic 
values, as being part of a mechanism for discrimination of detail. 

Let a light line segment on dark background be presented and 
let L be its geometrically similar image on the retina along which the 
photosensitive elements are excited. Let L, be a corresponding image 
on the central retina (Rashevsky, 1938, p. 266). Let the excited ele- 
ments along L, each send efferent processes in a random manner to 
a higher center (as discussed above) which is the group of synapses 
between the latter efferent ends and the afferent ends of fibers of the 
type discussed in section IT; and let these latter fibers be arranged in 
the order in which excitation occurs with increasing stimulus strength. 
If the afferent ends of each of these latter fibers (fiber pairs) are also 
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in synaptic contact with the efferent collaterals of a neural cireult 
showing periodicity then, under appropriate conditions, transmission 
will occur through the fiber pairs to a higher order fiber as though an 
adequate stimulus, of variable width, had been moved back and forth 
along the line of the afferent ends of the fiber pairs. If for a certain 
line L,, there is a specific linear center of the fiber pairs, we then 
have a possible mechanism exhibiting the properties suggested (Rash- 
evsky, 1938a). 

Summary. General fibers, satisfying the conditions discussed in 
section II, are used in connection with a mechanism which has the 
property of “following” a contour. 


IV. The reaction to the size of an object independent of its dis- 
tance. It is a well known phenomenon that the size of an object can 
be recognized within a considerable distance range (Kohler, 1929), 
and a great number of comparisons influence the judgment. However, 
let us assume that all cues are eliminated except those obtained from 
the image of the object on the retina and from the muscles of the 
eyes, the latter being in focus upon the object. N. Rashevsky (1938, 
p. 276) has discussed the possibility of interpreting the phenomenon 
under such conditions. Let us consider the neural mechanism of Figure 
4 which has the property of giving a representation of the size of an 
object irrespective of its distance. 

Let proprioceptive nerve fibers which arise at the fibers of a 
muscle of the eye be collected together into a single tract and let these 
fibers or their collaterals indirectly influence a group of neurons N. 
The nerve fibers will have some distribution of thresholds so that for 
various degrees of contraction or extension of the muscle, various 
numbers of fibers will be excited and the intensity of excitation, E , 
in each will vary with the degree of contraction. We shall hereafter 
refer to each tract of nerve fibers as though a single fiber and treat 
the group N as a single neuron. 

Let the proprioceptive nerve fibers I, IT, III, IV arise respectively 
from the lateral rectus, medial rectus, zonula ciliaris, and ciliaris mus- 
cles of the left eye and corresponding fibers (primed) to the right 
eye (Figure 4). Assume that the corresponding fibers from each eye, 
as I and I’, are equivalent. Let fibers J and III be of the simple in- 
hibitory type producing only 7 , and let fibers IJ and IV be of the sim- 
ple excitatory type producing only ¢. We shall also consider only the 
stationary state so that « or 7 is directly proportional to the total ex- 
citation, E’, in the fiber (Rashevsky, 1938, Ch. AALL): 

If the subscripts 1, 2, 3, 4 refer to the fibers J, IT, III, IV or the 
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FIGURE 4 


corresponding muscles, and if x; is the degree of contraction of mus- 
cle i from some position, then #;(%;) is the intensity of excitation in 
the i-th fiber when the corresponding muscle is contracted by an 
amount x;. Considering only those nerve endings which are excited 
by the contraction of the muscle, we have Hi(Yim) = 0, where Zim is 
the minimum degree of contraction producing excitation and is a 
threshold value, and also E;(x;) is some function which increases 
monotonically and is concave downwards. Equations (1) and (2) 
[Z'(S)] satisfy these conditions. Let us however assume that, over a 
sufficiently wide range, a linear relation holds in the case of fibers I 
and JI and that the fibers are equivalent except that one produces 7 
and the other ¢, OF 21m = Lom = Um and fb, = f. = B, and hence 


EB, (21) = B(%,— Cra )as (15) 
E,(%2) = B(%2— 4m) . (16) 


Let the eyes be focused upon an object which is “straight ahead” 
and let the distance from the eyes to the object of area s be d. Also 
let p be the distance from the fovea to the optical center of the lens 
system, 21 be the distance between the eyes, and 26 the angle, at the 
object, formed by the lines from the eyes. If o is the retinal size of 
the object, we have the relations 


o p? 
-=—=, 17 
BE fe a) 
and 
tan§=5 sind 0, (18) 


the latter approximations being accurate except for objects within a 


few inches of the eyes. 
Now for an object indefinitely far away let the states of contrac- 


tion be Z, = %» = &@ for the muscles corresponding to I and IT and let 
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Ax; be the change due to convergence of the eyes. Then we have 


= CALs (19) 
fo =¢-+ Aa, (20) 
and 
Ne (21) 
since 
x, + % = 4, -- % = constant, (22) 


or the degree of extension of one equals the degree of contraction of 
the other. But, since the change in degree of contraction is propor- 
tional to the angle and the latter is given by equation (18), 


At =—Am=kd=— (23) 


where k is a fixed parameter depending on the lengths of the muscles 
and the diameter of the eye. Substituting (23) into (19) and (20) 
and the results into (15) and (16), we have for each eye 


E, (t) = 6@—t%——=) » (24) 


E, (2) = BG — p+ Sy. (25) 


Now «, = A, E2/d. and 7, = B, E,/b, ; but as the fibers were 
assumed equivalent A./a, = B,/b,. Then at synapse s, we have for 
the next excitation (e — 7) from the fibers J, IJ, I’ and II’, of both 
eyes, 


AA, Bk  hsdy 
Ls Aegon tn Say 
where h; is the threshold of V , and dy), = 4 Az f kl/a.h;. Then let the 
fiber V be such that the intensity of excitation H; may be approximat- 


ed by a logarithmic function of the stimulus (Rashevsky, 1938, p. 264) 
so that we may write, if K; is a constant depending on fiber V , 


(s— 9) . d 
Ey = Ke log ——— = 5 log = . 


(e—J)1 = 2e,—2h,= (26) 


(27) 


As the eyes change accommodation from a distant to a closer ob- 
ject, the ciliaris muscle of each eye contracts, increasing the excitation 
E, in the excitatory fibers IV and IV’, and the zonula ciliaris of each 
eye is relaxed, decreasing the excitation EF in the inhibitory fibers III 
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and III’. The actual mechanism here is complex but the net excita- 
tion (e — 7). at synapse s, increases with the reciprocal of the dis- 
tance. Let us assume that the result is such that we may approximate 
the intensity of excitation E, in VJ by an equation similar to (27) or 


dos 
E, = K, log —*. (28) 


Let the contour of the retinal image of the object be represented 
at a center C, (Figure 4) as discussed by Rashevsky (1938, Ch. 
XXVII) so that the net excitation factor is proportional to the length 
of the contour. For similar retinal images the length is proportional 
to the square root of the retinal area, ¢. A particular value of the ex- 
citation at C; then represents a particular retinal size, the shape being 
constant. Now let an inhibitory fiber VII, from C; to N , also be of a 
type for which we may use the logarithmic relation between the ex- 
citation E, in the fiber and the net excitation factor, proportional to 


Vo , acting as the stimulus, so that if K, and o are constants, 


ee Kalo (29) 


do To 


E, = 2K, log 


Since e = A E/aand j = B E/b for the stationary state, we have 
from equations (27), (28), and (29), for the net excitation (e==9), 
= «,-+ e, — j, acting on N, 

5 K; daa A, Ke des B; K; 


f A o 
ial enti, <1 08 7 ae ‘a log 7 per e ae (30) 


If A,K;/d; = AcKe/@. = B,K,/b,, we then have for the above equa- 
tion, because of (17) 


A; K; dor doe Com A; K; dor doo do 


se og ee oe ee 10 31 
(e—J)n T log —y— ed (31) 
or 
: A; Ks hy 
as — —y> , oe 
(e—J)n = log — (32) 


where 8) = do doz o0/p? . Then under the above conditions, N receives 
stimulation which depend only on the variable, s, which is the size 
of the external object. 

It will also be noted that if the object is not directly in front of 
the subject but to one side so that both eyes are turned from being di- 
rected straight ahead, it can be seen that the net excitation from J and 
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II will still depend only upon the angle between the eyes if equations 
(24) and (25) are assumed. 

In order to arrive at equation (32) it was necessary to make a 
number of restrictive assumptions. Let us show that under much 
more general conditions we may obtain approximately the same re- 
sult. From the previous discussion it is clear that, on all counts, the 
values of H; and H, increase with the reciprocal of the distance d 
from the object to the eye. Let the total excitation factor acting upon 
N due to H; and E, be given by «,(d) which then increases with 1/d. 
Similarly the inhibitory factor 7, due to the center C,; increases in 
some manner with the linear retinal size but the latter, for a given 
object, varies directly with the linear size of the object and inversely 
with the distance d, and hence j;(d) is an increasing function of 


\/s/d. The net excitation acting on N is then 


(€ —J)n = &a—Jr.- (33) 

Let e, and 7, be given by the series 
coger ee (34) 
jy get, ee tea (35) 


where d, , d, and s) correspond to threshold values. Also let é and ” 
be defined by 


=s—s, (36) 


where s and d are values of s and d which are well removed from ex- 
treme values. Then substituting for s and d in terms of & and n in 
(34) and (35), we have upon expanding, 


&—jr=atpéEtynt--., (38) 


where 


: (39) 
= — 14 2 aera s : 
Z L 2+ as Bs] 
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In equation (38) only the first order terms in & and 7 from the 
expansions are written, and only those from the first three terms of 
(34) and (35). If we impose the restriction that 


VO. (40) 


we may then write for the net excitation acting on N, from (33), 
(36) and (38), 


(e¢—j)n =a4 B(s—s)+---, (41) 
so that, to a first approximation, the net excitation factor acting on N 
is a function of the object’s size, s, only. If an object of size s ata 
distance d is used as a standard, a similar object of size s and distance 
d, values not too removed from s and d, will influence N depending 
only on the relative sizes. The deviation from a judgment of size only 
will increase as either s or d deviate from the size or distance of the 
standard and will increase as the standard deviates from s and d. The 
size and distance of an object at which judgments are determined 
most completely by size alone are given by s and d. 
Summary. A neural mechanism has been discussed which has the 
property, under certain conditions, of giving a representation of the 
external size of an object independent of its distance. 


V. An interpretation of the effect of a warning stimulus and of 
stimulus strength upon reaction time. Let us set up a simple abstrac- 
tion of the situation in which a stimulus produces a fixed response but 
also in which another stimulus produces a “set”. One of the simplest 
possible abstractions is given in Figure 5 in which the stimulus S pro- 


FIGURE 5 


duces a response, R, and a warning stimulus S,, influences the re- 
sponse. 

Let the subscripts 1, 2, and 3 refer to the fibers J, II and 
III respectively. Also let E,(S) and E;(S,.) be the respective mono- 
tonically increasing functions of the stimuli, S > h, and Sy > hs. 
Since we are considering only values of S sufficient to produce the re- 
sponse without S,,, the net excitation (e — 7) due to I must increase 
monotonically from zero to a value h, at least. Hence without effect-' 
ing any essential change we may assume that I is a simple excitatory 
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fiber producing only ¢. In the case of fiber IJ we need not specify any- 
thing except that when the net excitation at synapse s exceeds hy. a 
response is produced after a time ¢) , a value which we shall assume to 
be essentially a constant and which includes the conduction time and 
the delay at the muscle. Fiber JJ we shall assume to be of the gen- 
eral type producing both e and 7. We then have, similar to equations 
(3) and (4), the following equations for the excitatory and inhibitory 
factors at synapse s (Figure 5), 


d &, 

TE Aine (42) 
d 

TE = 4sBs—t 60» (43) 
dj 

i = B, E,;—Dsjs> (44) 


where the nerves have different parameters, and where the times ¢, 
and f,; are not measured from the same origin. Since the E’s are con- 
stant we integrate and obtain the relations 


A, E 
& = = — (Leary, (45) 
A,;E 
eg = 2 (1 es), (46) 
as 
. BE 
js = (1— eh), (47) 


where the constants of integration are determined by the initial con- 
ditions e, = 0 for t, = 0 and e, = j, = 0 for t; = 0, so that the times 
t, and t, are measured from the times of presentation of S and S,, 
respectively, but do not include ¢, . 

The net stimulus acting on JJ is then 


6 Ag. 
1 cide Bld Seer - 6 gee) on 


(48) 


(1 ‘es e-ats) pa (1 pao e e-bsts), 


A; E; B; E, 
a; b 
The time t, between the presentation of the stimulus S and the initia- 
tion of the response, ¢’; -+ t, , is given by the solution of equation (48) 
for ¢, = @, when e,-+ ¢,—j,; = h., and adding t,. Solving for t’, , 

we have * 
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po a ay A, E 
cd pet wen Gans Mckas - 8 (1 — emt) 
(49) 
BE; bt! 
x (1 — ets) }]. 


In the above equation t’, is the time between the presentation of S, 
and the beginning of the reaction, less the time t,. If the preparatory 
interval, t, = t’; + t', + t., is much greater than the reaction time, 
as is generally the case, we may neglect t’, + t, and refer to @’; as the 
preparatory interval. However, for very short warning periods, the 
proper correction must be made. 

In order that S produce the response in a finite time and that S, 
alone should not do so, t’; must be positive and finite for all values of 
t’,. We then have the condition 


A,E 
hy > —— (1— e*") 1 Pe ie es) 
Gs Ds 
(50) 
pas A, E, 
ay 


Equation (49) gives the response time, t, = ¢'; ++ t, as a function 
of the warning period, t’,, the strength of the stimulus S in terms of 
E,(S), and the strength of the warning stimulus S. through E; (Sw), 
in terms of the parameters a,, A,, 43, A;,6;,Bs,h, hy, he, and ty. 
From the equation we find that the response time #, decreases with 
increasing stimulus intensity S such that for a minimum value of S 
there is no response and for large values of S the time t, decreases to 
a minimum value, since E,() is finite. Also t, decreases with the 
intensity of S,, from a finite initial value to a minimum value, if t’; is 
not too large, in which case ¢, may increase with S,, from its initial 
value and approach a limit. The effect of t's; upon the time t, is best 
illustrated by the curves below. It should be noted here, however, that 
the value of t, for t’; = 0 and that for ’; = o are not necessarily 
equal although they are so taken in the graphs. 

Data by H. Woodrow (1914) give the relation between t, and t's 
and the relation between ¢, and S. A rather complete check could be 
obtained but unfortunately the values of S were not measured. Hence 
separate sets of curves will be used to illustrate the above relations. 
; In order to determine t,(S) it is necessary to specify E,,(S). That 

the relationship is approximately logarithmic over a wide range for 
single fibers is indicated by data of H. K. Hartline and C. H. Graham 
(1932). Such a relation has also been used in previous theoretical 
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developments (Rashevsky, 1938), (Landahl, 1938). If we then intro- 
duce ‘in See 


E, = K log» S—K Ey, (51) 


where K E, takes care of the unit of intensity which is generally not 
given, we have for equation (49) 


1 H 
F eck waGaam, donteelond Ses 
or 
fi login S— Ey 
=t+— ——______- ye 
ty to = moe Fre EY 5 EET ] ? (5 ) 
where 
Cy A; E; 
—— h. ee a ee —Agt’s 
cf ALS Z a3 CLaiiede dia: 
(53) 
B, E, (1 = eh) } : 
b; 


Equation (52), which is similar to the first half of Rashevsky’s, 
equation (12) (Rashevsky, 1938, Ch. XXII), is compared with experi- 
ment as shown in Figures 6 and 7. The data in Figure 6 are by G. O. 
Berger (1886) and J. McK. Cattell (1886) for subjects B and C. The 


BERGER ANO CATTELL 
VISUAL DATA 


OSUSJECT 8 
@ SUBJECT C 


REACTION TIME IN SECONDS t, —= 


0 q 
| 2 1 B58 16 Fes: 100 
STIMULUS INTENSITY S—> 


FIGURE 6 
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stimulus intensities were measured photometrically. The range of in- 
tensity is from 1 to 1000 so that a logarithmic scale is used. The low- 
est intensity was near the threshold, being occasionally unnoticed. 
Each experimental point was the average of approximately one hun- 
dred and fifty reactions. A warning signal of a constant time interval 
was used. The units of intensity are not given, the smallest being 
given the value 1. For the upper curve the following values were 
assigned to the parameters in equation (52) to obtain the curve: 
to = 198 sec., 1/a, = .418'sec., Hy = —.85, H = .25 (subject B) ; and 
for the lower curve: f = .153 sec., 1/a, = .279 sec., Hy = —.68, 
H = .25 (subject C). 

The threshold, S’, is given by log S’ = E,+H. For subject B, 
the threshold assumed is S’=4, and for subject C, the threshold 
is S’ = 8. The actual thresholds are not given but are not much less 
than 1. The above values may possibly be too low. 

A large number of experiments under various conditions have 
been performed by H. Piéron (1920), most of which were visual re- 
actions. For most of the cases about twenty reactions were averaged 
for each value. Although the number of reactions is smaller for each 
value, the results are very significant in view of the wider range of 
intensities, in several cases, and the variety of situations investigated. 

Visual data by S. Froeberg (1907) when plotted on a logarithmic 
scale show a linear relationship. The range is from 1 to 100 and each 
value is an average of 400 reactions. The threshold was considerably 


PIERON: AUDITORY DATA 


REACTION TIME IN SECONDS | t, —» 


480) A ae 100 jeoe. 
STIMULUS INTENSITY S$ —= 


FIGURE 7 
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less than 1, the smallest value used. The linear relationship between 
the reaction time and the logarithm of the intensity is contradicted by 
the results of G. O. Berger (1886) and McK. Cattell (1886) and also 
by the numerous experiments of H. Piéron (1920), in one case cover- 
ing a range from 1 to 100,000. It should be noted that the curve must 
be concave upward for small stimulus intensities if a threshold exists 
or unless there is a discontinuity. Also, because of the method used 
to measure the intensity, a slight difference in background from abso- 
lute black would shift the measurements so as to make the resulting 
curve more like that in Figures 6 and 7. The large intensities would 
show no shift while near the smallest values used the shift would be 
of the order of 100 times the background intensity. 

In Figure 7 auditory reaction times are plotted against intensity 
on a logarithmic scale (Piéron, 1920). The curve is computed from 
equation (52) using the following values for the parameters: t) = .185 
sec., 1/a, = .304 sec., HE, = —.16, H = .10. The threshold is taken as 
S’ = .87, whereas unity is given as the threshold. 

In Figure 8 gustatory reaction times are plotted against inten- 


PIERON: GUSTATORY DATA 


STIMULUS INTENSITY S —= 
FIGURE 8 


sity (Piéron, 1920). The scale in this case is linear. Also, in order 
to fit the curve so as to include the first point, S = 1, it was necessary 
to assume a linear relation between E and S, that is, EF = B(S — h,) 
(Rashevsky, 1938, Ch. XXII, eq. 3), instead of equation (51) so that 
instead of equation (52) we have from (49), if H, = H K/s, 


Sh, 
Se (54) 


1 
| t, = ty + a, loge 
The curve in Figure 8 is obtained by the above equation where the 
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(lane have the values: t = .220 sec., 1/a, = 3.800 secs., hi, = .6, 
H, = .2. The value of S = 1 is given as the threshold ene 1920) 
Perens we have assumed a value of S’ = .80. 

The range for the gustatory data, Figure 8, is from 1 to 33. Un- 
fortunately the range could not be extended. The use of the linear 
relation between E and S then depends entirely upon the first point. 
Were it possible to extend the results for higher intensities and to 
insure the correctness of reaction times for low values of the inten- 
sity it would be possible to determine whether a logarithmic or linear 
relation between EF and S would be most accurate over a wide range. 
It might possibly be that all the relations are very nearly logarithmic 
or if not, only certain types are linear, as perhaps cases of direct 
chemical stimulation. 

Since the values of the reaction time for low intensities are so 
important the following suggestion may be considered. If, instead of 
the mean of a group of reaction times, we use the median, certain ad- 
vantages are obtained. For large stimuli and short times the mean 
and median will be very nearly equal. However, for small values of 
the intensity the situation will be different. For values even several 
times larger than the threshold no response may be made, or if made, 
may be excessively long. In either case the mean is not meaningful. 
If the median is used, it becomes possible to obtain time values almost 
up to the threshold value, the latter being then taken as a value for 
which the response is made in infinite time for just one half the trials. 
For any greater stimulus value the median is finite. Consider a re- 
sponse as having been made if half the trials or more are responded 
to or noticed and let the time be the median of the response times. 
Then the threshold is that value for which response is made in infinite 
time, and it corresponds to the definition in equation S’ = E,-+ H. 
On the other hand the mean is infinite unless all failures to react are 
arbitrarily eliminated. And one excessively high value entirely out- 
weighs the remaining values and it becomes difficult to define an ex- 
perimental quantity corresponding to the threshold in the equation. 

For a constant stimulus the t,(t’;) relation may be written from 


(49) 


bit —= log[M ++ J (ess — emt) ] (55) 
where 

= a (56) 

pa AsHats _ BBs (57) 
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It will be noted that in (57) we have assumed that t, is the same for 
t’, = 0 and t’; = 0, a restriction which is not necessary but which 
has some justification. 

In the experiments by H. Woodrow (1914) bearing on the above | 
relation three subjects were used. The data from one subject were 
chosen for illustration in Figure 9. The particular subject, V;, was 


WOODROW: SUBJECT Vs 


_ 


12 16 20 24 
PREPARATORY INTERVAL IN SECONDS ty = 


FIGURE 9 


chosen because the curve resulting from his performance is most typi- 
cal of the three (see Woodrow’s statement p. 24) and because the data 
were most complete. That the trend here is significant has been noted 
by other writers (Woodrow, 1914, p. 16). The two curves in the figure 
were obtained under conditions differing in one respect. For the case 
illustrated by the upper curve the subject did not know the length of 
the warning period whereas in the other case the subject was given 
practice with a particular warning period after which his responses 
with that warning period were recorded. 

The mechanism shown in Figure 5 has no hysteresial element and 
hence as such could not account for any learning. In fact, the warn- 
ing stimulus has to be assumed to be maintained throughout the inter- 
val. The addition of a simple conditioning circuit of either the single 
or double type (Rashevsky, 1938), (Landahl and Householder, 1938) 
would eliminate the restriction without necessarily complicating the 
results appreciably. However, the conditioning of the warning stimu- 
lus by the response is evidently of great importance. A more complete 
treatment of the problem would require a mechanism which, without 
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conditioning, would give the upper curve in the figure, and which, 
when conditioning is fairly complete, would give the lower curve. 

Pending such a solution let us complicate the mechanism just 
enough to justify the comparison of the theoretical properties of the 
mechanism with experiment under analogous conditions. Let an effer- 
ent collateral from III (Figure 5) lead to a conditioning center C. 
Let an efferent collateral from IJ also lead to C and from C let a fiber, 
IV, of the type of III lead to the synapse s so that S,, tends to be con- 
ditioned to the response. However, let the parameters be such that the 
net excitation at s due to both I/J and IV is insufficient to produce the 
response. Now if the warning periods, t,,, are of random length, no 
conditioning takes place and fiber IV carries no excitation. Equation 
(55) then may be used. If, however, a given warning period tends to 
be conditioned to the response, both III and IV influence synapse s. 
The additional excitation due to IV decreases the time reaction t,. If 
the effect of IV is much greater than III, the form of equation (55) 
may again be used but quite different parameters will hold. 

For the upper curve in Figure 9, the following values of the para- 
meters in equation (55) were used: ¢, = .130 sec., 1/a, = .200 sec., 
Mm = 638, J = 2.00, ad, = .132 sec*., b, = .120 sec*. For the lower 
curve in Figure 8 the values of the parameters in equation (55) were: 
fa — 2150.sec.,.1/a,. — <200. sece Mr— .670, J.= Al, Gy = 1:28 sec, 
b, = .128 sec?. The data were obtained by Woodrow (1914) from 
subject V,. Each of the upper set of points represents an average of 
100 reactions and each of the lower set represents an average of 150 
reactions. Note that t, and a, are taken equal for each case and that 
M and b; are nearly equal in each case. The value of 1/a, is of the 
same order of magnitude as that in the cases of visual and auditory 
stimuli, Figures 6 and 7. The upper curve could not be made to pass 
nearer the highest and last point. A change in a, , which would not 
be easily justified, might raise the curve at the right. However, the 
average curve for the three subjects does not show such a sudden rise 
between t,, = 16 and f, = .20. 

It is of interest to note that the value of 1/a, is of the same order 
of magnitude for the above curves (Figure 9) as that for the visual 
data (Figure 6) and auditory data (Figure 7), although there is no 
apparent connection between the curves. However, it should be noted - 
here that the intensity time curves are largely determined by the 
threshold S’ , the value ¢, and the product of H and 1/a,. As long as 
_H/a, remains constant changes in H or 1/a, do not appreciably affect 
the shape of the curve. Then making the values 1/a, the same for two 
different situations imposes a condition upon H. If then any one is 
determined independently, the others are fixed. 


118 MATHEMATICAL BIOPHYSICS 


Summary. The properties of a simple neural mechanism (Figure 
5) are developed theoretically and applied to an experimental situa- 
tion. The development leads to a relation [equation (49) ] among the 
experimental variables S , the stimulus intensity, ¢, , the reaction time, 
and t,,, the warning period. The relation between the reaction time 
and stimulus intensity, derived from the general relation, is compared 
with experimental data in Figures 6, 7 and 8. The relation between 
the reaction time and the warning period, also contained in the gene- 
ral relation, is compared with experiment in Figure 9. 

The author is indebted to Dr. A. S. Householder for checking the 
numerical and algebraic calculations and reading the manuscript; 
and to Professor G. von Bonin for making a correction in the first 
section. 
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